names(list_of_cents)
[1] "fold_change_aritmetic" "fold_change_geometric" "fold_change_quadratic" "fold_change_harmonic" "delta_aritmetic" "delta_geometric" "delta_quadratic" "delta_harmonic"
Density plot of centralities
plot_density <- function(df, titulo ){cut = 0.01
df %>% gather() -> df_gathered
df_gathered['value'] %>% filter(abs(value) > cut) %>%
ggplot( aes(x=value)) +
geom_density(fill="deepskyblue4", color="azure4", alpha=.8)+
theme(legend.position="top") +
ylab("Density") + ggtitle(titulo) +
xlab("Centrality") + geom_vline(xintercept=0, linetype="dashed", color = "red", size=.4) -> Density_plot_of_centralities
return(Density_plot_of_centralities)}
map2(list_of_cents, names(list_of_cents), plot_density)
$fold_change_aritmetic
$fold_change_geometric
$fold_change_quadratic
$fold_change_harmonic
$delta_aritmetic
$delta_geometric
$delta_quadratic
$delta_harmonic
#centralities <- Results[,!is.na(as.numeric(colnames(Results)))] %>% gather()
Density_plot_of_centralities.panel <- ggarrange( NULL,Density_plot_of_centralities, nrow = 2, heights= c(1, .4), labels = c('a','b'),
hjust = -0.1, vjust = .8)
Density_plot_of_centralities.panel
fba_list_of_cents
$fold_change_aritmetic
$fold_change_geometric
$fold_change_quadratic
$fold_change_harmonic
$delta_aritmetic
$delta_geometric
$delta_quadratic